Measuring order in the isotropic packing of elastic rods 
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The packing of elastic bodies has emerged as a paradigm for the study of macroscopic disordered 
systems. However, progress is hampered by the lack of controlled experiments. Here we consider a 
model experiment for the isotropic two-dimensional confinement of a rod by a central force. We seek 
to measure how ordered is a folded configuration and we identify two key quantities. A geometrical 
characterization is given by the number of superposed layers in the configuration. Using temporal 
modulations of the confining force, we probe the mechanical properties of the configuration and 
we define and measure its effective compressibility. These two quantities may be used to build a 
statistical framework for packed elastic systems. 

PACS numbers: 46.32.+X, 46.70.Hg, 61.43.-j 



Packed elastic objects are ubiquitous in Nature and 
technology. For instance, DNA is packed in the cell nu- 
cleous or in viral capsids [1, 2], while growing tissues can 
be confined by their environment [3|,13]- The optimiza- 
tion of folding is crucial in the design of self-deployable 
structures, such as tents or solar sails (Hf, or in waste 
disposal. Like a granular pile, a confined plate can be 
either in a crystalline state, the stacked facets obtained 
by repeatedly folding a sheet into two, or in a disordered 
state, exemplified by a crumpled ball [6]. When a sheet is 
confined, the number of metastable configurations blows 
up 0. Meanwhile, self-avoidance leads to jamming be- 
cause it prevents the system from exploring the space 
of configurations. This raises the question of whether a 
confined sheet can be viewed as a glassy system, in the 
same class as a static granular medium [8]. Indeed, theo- 
retical studies proposed thermodynamical approaches for 
packed rods [2|, |9|, llO|. On the experimental side, a diffi- 
culty in the study of crumpled balls [11-16] arises from 
the hand-generation of configurations. In this context, 
the confinement of a rod in a p lane was an important and 
useful simplification [7l, [l7l - [l9| . However, drawing general 
conclusions from these experiments can be questioned be- 
cause of issues such as friction between the periphery of 
the configuration and the container, the anisotropic in- 
jection of the rod in the container, plasticity or the im- 
possibility of unfolding a configuration. 

Here we reconsider the packing of a rod in a plane. 
We devised an experiment allowing us to reversibly con- 
fine a rod by a central force, deriving from an isotropic 
radial potential. As a consequence, there is no contact 
between the container and the periphery of the config- 
uration, while the intensity of the forcing is controlled 
through the stiffness of the potential. We investigate the 
emergence of geometrical order through the stacking of 
layers. As this setup enables the temporal modulation of 
the confinement, we also probe the mechanical properties 
of configurations, and define an effective compressibility 
of a configuration, which we associate with geometrical 



order. We thus obtain a coupled geometrical and me- 
chanical characterization of the system. 
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FIG. 1: The experiment, (a) A rod of density 1 is confined 
in a circular Hele-Shaw cell, filled with saturated salted water 
of density 1.16. The rod is in its initial configuration, before 
the cell has been set in rotation, (b) The cell is made of 
a superposition of three disks; the intermediate one (ii) is 
pierced with two holes to enable equilibration of pressure into 
the liquid. The chambers are made watertight with flat and 
toric joints shown in dark, (c) Example of a configuration 
obtained when a centripetal force is generated by the rotation 
of the cell around its axis, (d) Corresponding skeletonized 
image in which the number of layers per branch is defined. 

The principle of the experiment is as follows. A cir- 
cular Hele-Shaw cell is filled with a liquid and entrained 
by a motor. A rod is inserted into the cell. The cell 
is slightly thicker than the rod, so that the rod cannot 
cross itself, constraining a two-dimensional folding. The 
liquid is denser than the rod. As a consequence, when 
the cell is rotated, the rod is submitted to a centripetal 
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force and thus confined in a radial, parabolic potential. 
The cell was fixed in a vertical position on the axis of 
a 2kW motor, through a double ball bearing to avoid 
vibration transmission from the motor to the cell. The 
rotation velocity was controlled using an electronic fre- 
quency variator. We used a PDMS rod of circular sec- 
tion, = 2 ± 0.2 mm in diameter and L = 3 ± 0.01 m in 
length and density p = 1. The liquid was salt-saturated 
degassed water of density pi = 1.16 at ambient tempera- 
ture. The cell was initially made of two disks of 50 cm in 
diameter: a stiff one in lOmm-thick Duralumin (z), and 
a transparent one in 15 mm-thick Polycarbonate (n), al- 
lowing us to observe the confined rod (Fig. [TJi,b). Cham- 
ber 1 was filled with the liquid and the rod; the gap thick- 
ness was fixed at 2.5 mm by a ring of Plexiglas inserted 
between the two disks. However it turned out that this 
first design does not allow a constant gap. Indeed, the 
rotation induced a pressure gradient in the liquid, with 
a negative pressure at the center, leading to the bending 
of the Polycarbonate plate and the thinning of the gap. 
The solution was to add a second Polycarbonate disk 
{Hi) to obtain a second chamber of thickness 1.5 cm; this 
chamber is filled with the same liquid and communicates 
with chamber 1 with two holes of 1 mm in diameter, en- 
abling the equilibration of pressure in the two chambers. 
Chamber 2 serves as a sacrificial chamber: when the cell 
is rotated, disk (Hi) bends inwards (instead of disk (ii)) 
but the gap is thick enough so that it does not close at 
our maximal rotation velocity of 20 Hz. The chambers 
were filled through holes pierced in the edge of disks {i) 
and (iii) and the rod was gently inserted using a hole in 
the back of disk {i). The whole setup was placed in a dark 
room and lit with three stroboscopic lamps with diffus- 
ing screens. Movies were taken with a CCD camera. The 
duration of a fiash was short enough (1 ms) to get sharp 
images of the rod, even at the higher velocities. Using 
a computer interface, the camera and stroboscopes were 
triggered with the same square periodic signal, while the 
variator was controlled using a DC voltage. To enhance 
contrast, we used a white rod on a dark back: an adhesive 
sheet of black plastic was laid on the Duralumin disk (i). 
Binary images of the rod were obtained by thresholding. 

The initial configuration of the rod (Fig. [TJi) is pre- 
pared using 8 magnetic beads of diameter 1.8 mm in- 
serted into the first chamber and moved from outside 
with a magnet. The cell is then set in rotation; the 
time to reach the desired frequency is approximately 3 s. 
The control parameter is the frequency / of rotation (See 
EPAPS Document No. [number will be inserted by pub- 
lisher] for a movie of a typical experiment.) First, con- 
figurations can be characterized using the radius of gy- 
ration. 




where s is the curvilinear coordinate of the rod and r{s) 



the distance to the cell axis. This quantity is directly cal- 
culated using the binary image of the confined rod. The 
radius of gyration decreases with time, rapidly when the 
rotation is started (Fig. [2^) and then more slowly, reach- 
ing a plateau value in a time lapse of the order of 10^ s. 
This final value is not unique and differs according to re- 
alizations at a given frequency. As we only investigated 
equilibrium configurations, we had to wait 1800 s for each 
value of the frequency before taking measurements, to 
ensure that equilibrium is reached. When averaging over 
realizations, the radius of gyration is found to be a de- 
creasing function of confinement strength (Fig. [JJ)). It 
appears that, when the same experiment is repeated, a 
large diversity of sizes and geometries is observed. The 
aim of the present study is to quantify geometrical order 
in equilibrium states. To do so, we extract the skeleton of 
a folded configuration from binary images. Vertices are 
detected as self-contact points, ie points of the skeleton 
having three neighbors (Fig. [Th, d). We define branches 
as portions delimited by two vertices. A given branch 
may contain several layers of the rod: the thickness of 
the branch on the binary image directly yields the num- 
ber of layers. 
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FIG. 2: The radius of gyration, (a) The radius of gyration Rg 
as a function of time for three realizations in which the cell 
was launched from to 14 Hz in 3 s. Radii reach a plateau 
value in about 10^ s. The final radius differs according to 
realizations, (b) On average, the radius of gyration decreases 
with the strength of confinement (quantified by the rotation 
velocity, / of the disk). The error bars are sized according to 
the standard deviation of measured radii for each frequency. 

Inspired by observations (Fig. [3^), we first consider a 
geometrical definition of order using the number of su- 
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FIG. 3: Stacking, (a) Order and disorder. Two extreme ex- 
amples of equilibrium configurations obtained at two different 
frequencies, with a large and a small number of superposed 
layers, respectively, (b) Average number of layers per branch, 
Ni/ij, as a function of the radius of gyration Rg. Each point 
corresponds to one realization of the experiment. The colors 
of the symbols correspond to the imposed rotation velocity of 
the disk. An ordered configuration (large Ni/i) has a small 
radius of gyration. 



perposed layers. Indeed, the configuration of absolute 
minimum of energy is a spiral [7] , which can be qualified 
as very ordered and in which all layers are superposed and 
contained in a unique branch. Experimentally, this num- 
ber of superposed layers is defined as the average number 
of layers per branch in a given configuration. We jointly 
measured the average number of layers per branch TV^/^ 
and the radius of gyration of a configuration, as shown 
in Fig. [3]3. The data roughly collapse on a single curve, 
independently of the confinement strength. Therefore, 
either the radius of gyration or the number of layers ap- 
pear as better characterizations of the configuration. A 
disordered configuration, ie with small has a larger 
radius than a spiral ordered configuration. As in some 
previous experiments 



dis- 
tinctive feature of the confinement of rods, and our setup 
allows us to show that stacking decreases with the radius 
of the configuration. In the following we investigate the 
possibility of characterizing geometrical order without a 
detailed knowledge of the geometry of configuration. In 
other words, we seek an independent measure of order. 

We realized annealing experiments in which the con- 
finement strength was repeatedly increased then de- 
creased by varying the rotation velocity appropriately. 
At each step in rotation frequency, we waited 1800 s in 
order to reach equilibrium (Fig. [2^), so that an anneal- 
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FIG. 4: Compressibility, (a) Evolution of the radius of gy- 
ration Rg during an annealing experiment. After a first ir- 
reversible brach, the radius follows a reversible branch (par- 
allel to the dashed line). A configuration is characterized by 
the slope of the line, cind its intercept, Rq, with the axis 
/ = 0. (b) The characteristic radius Ro as a function of com- 
pressibility X- Each point corresponds to one realization of 
an annealing experiment. The error bars correspond to the 
uncertainties estimates of the linear fits as described in (a). 



ing experiment typically took 12 hours. After a few fre- 
quency steps, the radius of gyration followed approxi- 
mately the same line (Fig. |4^). The first steps are irre- 
versible, while the line is reversible. In other words, the 
system follows an irreversible branch in the (f^Rg) phase 
space before falling on a reversible branch. This behavior 
is reminiscent of the evolution of the density of a tapped 
granular pile according to the tapping acceleration, as 
reported in [20^. In our case, each annealing experiment 
can be characterized with the intercept, Rq, and with 
the slope, x, of the reversible branch. Furthermore, we 
observed that no topological changes occurred along the 
reversible branch, as illustrated in Fig. [5l the relative 
positions of loops remain constant, and the configuration 
only seems to breath. As a consequence a given annealing 
experiment leads to a well-defined configuration, which 
can be characterized with Rq and x- corresponds to 
its effective radius, while x defines an effective compress- 
ibility of the configuration as it measures its susceptibil- 
ity to the confinement strength. When plotting these two 
quantities, the characteristic radius Rq is found to be an 
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increasing, roughly linear function of the compressibility 
X (Fig. HJd). The results are independent of the confine- 
ment (frequency) since, whatever annealing path, points 
seem to collapse on the same master curve. As we found 
above that an ordered configuration has a small radius 
of gyration (Fig. [Sb), this second set of results indicates 
that an ordered configuration has a small compressibility, 
and, conversely, a disordered configuration is highly com- 
pressible. A possible interpretation is that self- avoidance 
imposes a stringent constraint on ordered configurations 
for which more stacking means less freedom in exploring 
phase space; as a consequence an ordered configuration 
would be less compressible. 




FIG. 5: Example of the evolution of a configuration during 
an annealing experiment. Rotation frequency was varied from 
12 Hz (left) to 6 Hz (right), by steps of 2 Hz of 1800 s duration. 
Each colored area corresponds to a loop. The relative position 
of loops is invariant. The configuration only breathes under 
the perturbation. 

To summarize, we built an experiment allowing the 
two-dimensional confinement of a rod in a parabolic po- 
tential. This experiment allowed us to quantify order 
and disorder in a configuration, using a geometrical mea- 
sure of stacking and a mechanical measure of (effective) 
compressibility. We found that these two quantities are 
strongly correlated with the characteristic radius of the 
configuration. Although this effective compressibility is 
strictly a susceptibility to changes in rotation velocity, 
it can be readily generalized to other systems such as 
crumpled balls. In that context, a salient feature of com- 
pressibility is that its measurement does not require a full 
knowledge of the geometry of the configuration. In other 
words, order could be inferred from the effective stiff- 
ness of the crumpled ball. Future work should address 
whether these quantities could be used in a thermody- 
namic approach to packing, or how these quantities could 
emerge in such an approach. Thus geometrical and me- 
chanical properties appear as strongly entangled in the 



packing of sheets and rods. 
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